Spontaneous population imbalance in two-component Bose and Fermi gases 
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We study two-component (or pseudo-spin-^) Bose or Fermi gases in one dimension, in which 
particles are convertible between the components. Through bosonization and numerical analyses 
of a simple lattice model, we demonstrate that, in such gases, a strong intercomponent repulsion 
induces spontaneous population imbalance between the components, namely, the ferromagnetism 
of the pseudo spins. The imbalanced phase contains gapless charge excitations characterized as a 
Tomonaga-Luttinger liquid and gapped spin excitations. We uncover a crucial effect of the intercom- 
ponent particle hopping on the transition to the imbalanced phase. In the absence of this hopping, 
the transition is of first order. At the transition point, the energy spectrum reveals certain degen- 
eracy indicative of an emergent SU (2) symmetry. With an infinitesimal intercomponent hopping, 
the transition becomes of Ising type. We determine the phase diagram of the model accurately and 
test the reliability of the weak-coupling bosonization formalism. 

PACS numbers: 05.30.Jp, 05.30.Fk, 03.75.Mn, 75.10.Jm 



I. INTRODUCTION 

Ultracold atomic gases offer highly controllable labo- 
ratories for testing and exploring novel many-body phe- 
nomena in interacting systems [[y, . One line of current 
interest is to confine atoms to highly elongated traps, ef- 
fectively creating one-dimensional (ID) systems. In ID 
interacting systems, the elementary excitations are col- 
lective modes, and the intuition based on the free-particle 
theory breaks down. For the spinless (one-component) 
case, theory predicts the equivalence of Bose and Fermi 
gases; both are described by the Tomonaga-Luttinger liq- 
uid (TLL) theory at low energies ||. As a hallmark ex- 
ample, fermionization of bosons has been observed in a 
Bose gas of 87 Rb atoms tuned into a strongly repulsive 
regime j|, ||. Another frontier of activity is the creation 
of multicomponent gases using different internal states 
of atoms or using different species of atoms. If we limit 
our attention to ID systems, a two-component Bose gas 
composed of the two hyperfine states of 87 Rb has been 
confined in a ID trap || Even without internal states, 
one can load the atoms in a double- channel trap j? 10 
[Fig. |l|(b)] or on a ladder-type lattice |ll], |Lq |, effectively 
creating a ID two-component gas. These two-component 
gases are expected to display a variety of phases depend- 
ing on the intra- and intercomponent interactions. In 
a double-channel or ladder structure, in particular, the 
magnitudes of these two interactions would be different 
essentially and controllable separately. 

In this paper, we study the two-component Bose or 
Fermi gases in one dimension and analyze a quantum 
phase transition induced by a strong intercomponent re- 
pulsion. In two-component gases consisting of two species 
of atoms, it has been argued that a strong intercompo- 
nent repulsion induces the phase separation (the demix- 
ing) of the species (see, e.g., Refs. |l3|— . Here we 
consider an analogous instability in the different situa- 
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FIG. 1: (Color online) Illustrations of ID two-component 
gases. The two components represent two internal states of 
atoms [shown by different colors in (a) and (c)]. Even without 
internal states, a double-channel trap potential can produce 
a similar situation [(b) and (d)]. A strong repulsion between 
the components induces the population imbalance [(c) and 



tion where the two components represent two (internal) 
states of single-species atoms and thus particles are mu- 
tually convertible between the components. The system 
has a Z2 symmetry with respect to the interchange of 
the two components. In this case, it is expected that, 
under a strong intercomponent repulsion, a single com- 
ponent can dominate the whole system (population im- 
balance) as in Figs. 0(c) and |](d). The Z2 symmetry is 
spontaneously broken while the translational symmetry 
is retained. This ordering may also be viewed as ferro- 
magnetism if we identify the two components with the 
pseudo spin |. A notable point in the present setting 
is that differently from the case of two-species mixtures, 
an intercomponent particle hopping can exist and vio- 
late the separate conservation of particle number in each 
component. It is also worth noting that the mechanism 
of spontaneous imbalance is in fact ubiquitous in nature. 
It underlies the ferromagnetism in [/(l)-symmetric itin- 
erant electrons [EOf , the vector chiral order in some frus- 
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trated magnets [^T], 22 1, and the spontaneous rotation of 
a trapped Bose gas [23|. This mechanism can occur both 
on a lattice and in a continuum since it does not involve 
any crystallization. 



II. MODEL 

We start from a two-component Bose gas model on a 
ID lattice defined by the Hamiltonian 



In spite of their simplicity and ubiquity, the basic 
properties of the spontaneous imbalance and related phe- 
nomena have been rather poorly understood. Mean-field 
analyses of coupled Gross-Pitaevskii equations can- 
not capture the effects of strong quantum fluctuations in 
one dimension. A study beyond the mean-field argument 
has been done based on the weak-coupling bosonization 
formalism jl4| . In this formalism, each component is de- 
scribed as a TLL, and the intercomponent coupling is 
perturbatively treated. Then, as will be explained in 
Sec. Ill, the spontaneous imbalance is predicted to occur 



at a point where one of the TLL parameters diverges and 
this formalism breaks down JTJ]. Therefore, this formal- 
ism cannot be used to describe the imbalanced phase nor 
the transition to it. 



Recently, ID spin-polarized (ferromagnetic) Bose gases 
in a strongly repulsive regime have been studied actively 
using integrable models [£4[ |25| and effective field theo- 
ries |26|-|29f . In these studies, the original Hamiltonian 
has SU(2) symmetry in terms of the (pseudo) spins, and 
this symmetry is spontaneously broken. In this setting, 
the excitations consist of a gapless charge mode char- 
acterized as a TLL and a gapless spin wave mode with 
a quadratic dispersion. In contrast, in our present set- 
ting, the (pseudo) spin SU(2) symmetry is reduced to Z2. 
The low-energy excitation structure in such a reduced- 
symmetry case has not been addressed. Furthermore, 
these studies of S'[/(2)-symmetric systems focus only on 
the properties of the fully polarized state, and the na- 
ture of the ferromagnetic (population-imbalance) transi- 
tion has not been discussed. 



In this paper, we analyze a simple lattice model of a 
ID two-component Bose or Fermi gas to address the basic 
properties of the spontaneous population imbalance. We 
use two numerical methods, exact diagonalization and 
infinite time-evolving block decimation p(| , in efficient 
manners to go beyond the existing effective theories. We 
determine the phase diagram accurately and test the re- 
liability of the bosonization prediction. Using the scaling 
of the entanglement entropy, we demonstrate that the 
low energy physics of the imbalanced phase is described 
by a one-component TLL, indicating that the excitations 
consist of a gapless charge mode and a gapped spin mode. 
We uncover a crucial effect of the intercomponent parti- 
cle hopping on the nature of the transition. Perturbation 
theory from the strong-coupling limit in the half-filled 
case provides a qualitative understanding of the phase 
diagram and the nature of the transition. 
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where 6 rj is a bosonic annihilation operator at site j in 
the rth component, and n r j — b], jb r j is the number 
operator defined from it. The first and second lines rep- 
resent hopping and potential terms in each component, 
and the third represents those between the components. 
We set t > and t' > 0. (The choices of the signs are ar- 
bitrary under gauge transformation.) For simplicity, we 
consider the hard-core limit Uq — > 00, where two parti- 
cles in the same component cannot occupy the same site 
j (but those in different components can) . 

We are also interested in the two-component fcrmionic 
model which is defined by replacing all the bosonic op- 
erators b r .j in Eq. ([!]) by fermionic operators f r .j. In the 
fermionic model, the hard-core interaction Uq automat- 
ically drops out. When t' = 0, the hard-core bosonic 
model is equivalent to the fermionic one via the Jordan- 
Wigncr transformation: 
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where the "string" part runs first in the first component 
and then in the second component. Because of this cor- 
respondence, the bosonic and fermionic models can be 
analyzed in parallel for t' — 0. At the special point 
t' = V = 0, the model (|l|) is equivalent to the solvable 
fermionic Hubbard chain, where the labels r — 1 and 2 
are identified with spin-up and spin-down states, respec- 
tively. In this case, the population imbalance is known 
not to occur. 

When t' ^ 0, the bosonic and fermionic models are 
no longer equivalent. For the following reasons, however, 
they are expected to display essentially the same physics. 
Specifically, we consider a different Jordan- Wigner trans- 
formation: 
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where the "string" part now runs alternately between the 
two components. Under this transformation, the intra- 
component hopping terms are transformed as 

bljhj+t = e*™-*fy id+1 , (4a) 
^>2 J+ i=e ! ™^+7 2 V 2j+ i, (4b) 

while other terms in Eq. ([j]) retain the same form. In the 
low- (respectively high-) density limit, the phase factors 
e '7r«r,! m gq_ (Q) are fixed essentially at unity (respec- 
tively — 1). At half- filling and for strong intercompo- 
nent repulsion, the intracomponent hopping t does not 
contribute in the first-order perturbation theory, as ex- 
plained in Sec. |v|. Therefore, at least in these cases, the 
bosonic and fermionic models lead to the same physics. 



III. WEAK-COUPLING THEORY 

We formulate a weak-coupling bosonization theory for 
the hard-core boson model ([j]) , and discuss the instability 
of the two-component TLLs following Ref. [l4| When 
t' = U = 0, the model decouples into two independent 
Bose gases, each equivalent to a solvable spin-i XXZ 
chain in a magnetic field. For —2 < V and < (n r ,j) < 1, 



each component r(= 
Hamiltonian B 



1,2) obeys a TLL described by the 



Hf = J Ax 1 - [K{d x 9 r f +K- 1 (d x <p r f] 
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where 9 r and <f> r are a dual pair of scalar fields, and 
x = J a o with ao being the lattice spacing. The group 
velocity v and the TLL parameter K can be determined 
from Bethe ansatz ||. To treat tf and U terms as per- 
turbations, we use the following bosonization formulas: 

n r ,j =Po~ a d x <p r (x)/y/Tr-\ , (6a) 

bl 3 = exp[-iV^r(x)](Bo + •••), (6b) 

where po = (n r _j) is the averaged density and -Bo is a 
nonuniversal constant. Introducing new bosonic fields 
<f>± = (fa ± </> 2 )/v / 2 and 0± = (6>i ± 9 2 )/V2, we obtain 
the effective Hamiltonian for Eq. (|l|), 
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with 



\ ttv J \ nv 

(8) 

Here we notice the effective separation of the two sectors, 
(<j)+,d+) and (</>_,#_), which are respectively called the 



"charge" and "spin" sectors by analogy with the Hubbard 
chain. In the weak-coupling regime (0 < P < v/ao), 
both the sectors remain gapless for t' = 0. Finite t' ^ 
opens a gap in the spin sector since the vertex term 
cos(v / 2~7r6U) with scaling dimension 1/(2KS) is always 
relevant. As U increases, these estimates (g) indicate 
v_ —> and K_ — > oo at 

U c = nv/(Ka ). (9) 

In other words, the coefficient V-KZ 1 of [d x (j)-) 2 in 
Eq. (Q) changes sign at this point. This indicates the 
breakdown of the bosonization description in the spin sec- 
tor. If we naively assume the existence of a term (9 2 .(/)_) 4 
with a positive coefficient in the effective Hamiltonian, a 
population-imbalanced state with An = (nij) — (n2,j) ~ 
—ao\/2/Tr(d x (f>-) ^ is expected to appear along with 
the breakdown po| , |23j . Equation (^ ) gives a naive esti- 
mation of the transition point [[TJ], |d ] . In the case where 
a two-component gas consists of two-species atoms, this 
breakdown corresponds to the demixing instability fl4| . 



IV. NUMERICAL ANALYSES 

To go beyond the weak-coupling theory and to address 
the strong-coupling regime, we employed two numerical 
methods for the bosonic model (Q): exact diagonaliza- 
tion (ED) and infinite time-evolving block decimation 
(iTEBD) |Fj. The iTEBD method generates the ground 
state of an infinite system through the use of a matrix 
product state. The precision improves as we increase 
the matrix dimension \- To perform a calculation at 
a fixed filling, the chemical potential /i was iteratively 
tuned through a feedback control in each iTEBD step. 
To achieve better convergence of the order parameter, 
we first performed the iTEBD for large U , where the or- 
der parameter is large. Then we gradually decreased U 
repeatedly using the obtained state as the initial state for 
the next U. We set t = 1 as the energy unit hereafter. 



A. t' — case 

Let us first analyze the case with zero intercomponent 
hopping t' = 0, in which the bosonic model (|l|) is equiv- 
alent to the fermionic one with the same Hamiltonian. 
In this case, the particle number of each component, 
N r = Ylj n r,jj is a good quantum number. Therefore, 
for a given total particle number N — N\ + N2, ED 
can be performed separately for each AN = Ni — N2 = 
0, ±1, . . . , ±N. The lowest eigenenergy in each sector is 
plotted in Fig. ||(a). We observe a direct change of the 
ground state from a uniform state AiV = to a fully 
imbalanced state AiV = ±iV as we increase U. Such an 
abrupt change of AiV indicates a first-order transition. 
All the energy levels cross at the same point U c = 9.34, 
which suggests an emergent SU(2) symmetry at the tran- 
sition point. Such an emergent SU(2) symmetry has 
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FIG. 2: (Color online) Numerical results for t' = 0, t = 1, 
and fHling=0.4. (a) The lowest energy levels for different AN 
obtained from ED and the density difference An from iTEBD 
(with x = 50), as functions of U. In ED calculations, the finite 
cluster of chain length L = 10 was used, (b) An obtained from 
iTEBD as a function of V. 



also been predicted in renormalization group analysis 
p9[. The recent strong-coupling theories for the SU(2)- 



9 are expected to apply at 
E(b), the density difference 



symmetric case in Refs. |26|-[2 
this point. In Figs. ||(a) and 
An evaluated by iTEBD shows a jump at a certain point, 
which also indicates a first-order transition. The transi- 
tion points obtained from ED and iTEBD are slightly 
different, which could be mainly attributed to an inher- 
ent hysteresis in iTEBD around a first-order transition 
point. The level-crossing point in ED shows only a very 
small dependence on the system size except when the to- 
tal density is close to zero or unity. Therefore, ED gives 
the better estimate of the transition point. 
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FIG. 3: (Color online) Numerical results for t' = 0.05, t = 1, 
and filling=0.4. (a) Density difference An and half-chain en- 
tanglement entropy Shaif obtained from iTEBD, as a function 
of U. The onset of An and the peak of 5*haif occur almost si- 
multaneously, (b) (An) 8 versus U, and (c) (An) 8 versus V 
near the transition point. The black solid lines in (b) and (c) 
are the linear fits of the data in the imbalanced phase. 



B. i! / case 

Now we analyze the effect of intercomponent hopping 
t' 7^ 0, in which Ni 2 are no longer conserved separately. 
As shown in Fig. |3|(a), the order parameter An calcu- 
lated with iTEBD grows continuously as a function of U, 
indicating a second-order transition. The peak of half- 



chain entanglement entropy Shaif [see Eq. (5) of Ref. |34| 
for its definition] gives a reasonable estimate of the transi- 
tion point. This quantity is known to diverge at a critical 
point H^, p3| , although finite \ introduces a cutoff to the 
divergence[B4| . In Figs. |(b) and |(c), (An) 8 is plotted 
as a function of U and V. The data are well fitted by 
a linear function except in the very close vicinity of the 
transition point where x dependence occurs. This result 
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indicates the relation An oc (17 - C/ c ) 1/8 [(V c - V) 1 '*] 
along the U [V] axis, in agreement with the the critical 
exponent /3 = 1/8 in the 2D Ising universality class |33] . 



C. Phase diagram 

The ground-state phase diagram is summarized in 
Fig. ||. The left and right sides of the phase boundary 
correspond to the uniform TLL and imbalanced phases, 
respectively. The transition points U c were determined 
by using the level-crossing point for t' = and the peak 
°f SWf for t' ^ 0. In Fig. |I](a), the phase diagram is 
symmetric under n — > 1 — n because of the particle-hole 
symmetry in the hard-core model. It is found that U c 
is shifted to larger values with increasing t 1 . Namely, 
the intercomponent hopping diminishes the imbalanced 
state. Figure |](b) shows the phase diagram in U-V 
space. We observe U c — > as V — >• — 2. This is natu- 
rally expected since, in the decoupled case (f = U = 0), 
the point V/t = — 2 corresponds to the two ferromag- 
netic Heisenberg chains in the spin-system language with 
S*j = n r j — l/2. The bosonization prediction (|9|) and nu- 
merical data agree well in this limit. On the other hand, 
as V is taken to zero, U c deviates from the bosonization 
prediction and tends to diverge. For V > 0, we did not 
observe a population imbalance, although the bosoniza- 
tion prediction (||) still indicates its occurrence. (We 
again note that the occurrence of the population imbal- 
ance can be disproved in the integrable case V = t' = 0.) 
This indicates intricate roles of U and V on the change 
of the TLL parameter A"_ , which are not covered in the 
weak-coupling approach. 
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FIG. 4: (Color online) (a) Phase diagram in the [/-filling 
space for V — —0.5. We denote the phase boundaries by 
solid, dashed, dotted, and dash-dotted curves for t' = 0, 0.02, 
0.05, and 0.1, respectively, (b) Phase diagram in U-V space 
for fillings 0.4. The Bosonization prediction (^) of the phase 
boundary is also shown for comparison. 



PERTURBATION THEORY IN THE 
STRONG-COUPLING REGIME 



where T z and Tf are pseudo-spin- \ operators defined as 



To gain a deeper understanding of the strong-coupling 
regime, we formulate a perturbation theory in the half- 
filled case (nij) + (n2,j) = 1. In this special case, the 
charge sector is gapped out and we can thus single out 
a simple structure in the spin sector. (This procedure 
is analogous to the derivation of the Heisenberg model 
from the Hubbard model for a strong on-site repulsion.) 
In the limit t/U,t' /U,V/U —> 0, the system decouples 
into independent sites. Each site j has doubly degen- 
erate ground states, | = \nij = l,n 2j = 0) and 
| JJ.)j = |m,j = 0, 7i2,j = 1). We use these as the basis 
of the Hilbert space, and we treat t, t' , and V terms as 
perturbations. First-order perturbation theory is equal 
to acting the projection operator P = Ylj(Hf}jj{1t\ + 
) j j (-U- 1) on both sides of t, t', and V terms in Eq. ([l]). For 
both the bosonic and fermionic models, the same effective 
Hamiltonian is obtained as 



(i) 



3 



(2VT*T Z +1 - 2t'T; 



(10) 



(ii) 



For t' = 0, Eq. © is a classical Ising model. A 
first-order transition at V — separates ferromagnetic 
(V < 0) and antiferromagnetic (V > 0) phases. The for- 
mer is nothing but the population-imbalanced phase. For 
t' > 0, Eq. (M) is equal to an Ising model in a transverse 
field, which is still solvable. Second-order transitions 
of Ising type separate ferromagnetic (V < —2t'), disor- 
dered (— 2t' < V < 2t'), and antiferromagnetic (2t' < V) 
phases. The presence of t' abruptly changes the nature 
of the transition, consistent with the numerical results 
in Figs. ^| and |^. Furthermore, the ferromagnetic phase 
diminishes as we increase t' , in accordance with Fig. |J. 
The appearance of the antiferromagnetic phase is due to 
a lattice effect specific to the half-filled case, and here we 
do not discuss it further. 
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FIG. 5: (Color online) Finite-interval entanglement entropy 
S L calculated by iTEBD. We set t' = 0.05, V = -0.5, and 
filling=0.4. (a) In the imbalanced phase, we obtain c ~ 1. (b) 
Approximately at the transition point, we obtain c « 1 + 1/2. 



VI. LOW-ENERGY PROPERTIES OF THE 
IMBALANCED PHASE 



We analyze the low-energy properties of the imbal- 
anced phase. In the case of t' = 0, it is clear that the low- 
energy physics is governed by a TLL in the charge sector 
since the fully polarized state is realized (see Fig. ||). On 
the other hand, the case of t' y= deserves to be inves- 
tigated. In the latter case, it is natural to consider the 
following two possibilities: (i) The effective separation of 
the charge and spin sectors still holds; the former is de- 
scribed by a TLL while the latter is gapped, (ii) The first 
chain (r = 1) with a dense particle density and the second 
one (r = 2) with a thin density separately form TLLs. 
To see which possibility is correct, we examine the scal- 
ing of the entanglement entropy Sl of the ground state 
|\&). For a finite interval of length L, it is defined as 
S L = -Trpn lnpo, where pa = Ttq|^)(*| is the reduced 
density matrix on Q obtained by tracing out the exterior 
O. In ID critical systems, this quantity enables one to 
determine the central charge c (an indicator of the num- 
ber of gapless modes) through the formula of a universal 



S L n~lnL 



const. 



(12) 



Figure |^ shows Sl calculated by combining iTEBD with 
the transfer matrix technique. Here, an initial state with 
a fixed number of particles was adopted to achieve bet- 
ter convergence of the wave function. In the imbalanced 
phase [Fig. ||(a)], c~ lis obtained by fitting the numeri- 
cal data with the scaling formula ( |i"2| ) . This supports pos- 
sibility (i). Namely, the dense and thin components co- 
operatively form a single TLL. This is strikingly different 
from a population-imbalanced state extrinsically derived 
by simply applying a magnetic field (chemical potential 
difference) for pseudo spins Figure ||(b) shows the 
result approximately at the second-order transition point 
U c for t' 7= 0. Although the data show some x depen- 
dence for large L, the fitting using almost ^-independent 
data for small L yields c ~ 1 + 1/2. This result indi- 
cates that an Ising-type transition with c = 1/2 occurs 
in the spin sector while the gapless mode with c = 1 in 
the charge sector remains intact. 



VII. CONCLUSIONS 

We have studied a simple lattice model ([[]) of a two- 
component quantum gas in the strong-coupling regime, 
using numerical and analytical methods. We demon- 
strated that the spontaneous population imbalance (fer- 
romagnetism) occurs as we increase the intercomponent 
repulsion U. We completed the accurate phase diagrams 
in Fig. |]. It was found that the spontaneous imbalance 
occurs only for V < and that the intercomponent hop- 
ping t' diminishes the imbalanced phase. While these 
intricate roles of V and t' cannot be covered in weak- 
coupling bosonization theory, perturbation theory from 
the strong-coupling limit for the half-filled case provides 
simple qualitative explanations of them. 

We have also uncovered the basic properties of the im- 
balanced phase and the transition to it. Using the scal- 
ing of the entanglement entropy, we demonstrated that 
the low-energy property in the imbalanced phase is gov- 
erned by a one-component TLL, indicating the separa- 
tion of gapless charge and gapped spin modes. This is in 
sharp contrast to the integrable fermionic Hubbard chain, 
where both the charge and spin sectors behave as TLLs 
even in the strong-coupling regime p7| . The transition 
to the imbalanced phase is of first order when t' = 
and of Ising type when t' 7= 0. In spite of this drastic 
effect in the spin sector, the t' term does not spoil the 
gapless property of the charge sector. At the first-order 
transition point for t = 0, the energy spectrum reveals a 
certain degeneracy indicative of an emergent SU{2) sym- 
metry. 

Although our analyses are done mainly for the bosonic 
model, we expect that the fermionic model also displays 
essentially the same physics as we discussed in Sec 0. 
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For t' = 0, the bosonic and fermionic models are exactly 
equivalent. For t' ^ 0, the t' term in the fermionic case 
has a different bosonization expression and thus a differ- 
ent scaling dimension from the bosonic case in Eq. (Q). 
We expect that, when this term is irrelevant (respectively 
relevant), the transition between the uniform TLL and 
imbalanced phases is of first order (respectively of Ising 
type). In particular, even in the presence of the t' term, 
the bosonic and fermionic cases become asymptotically 
equivalent (i) in the low- and high-density limits and (ii) 
in the strong-coupling limit in the half-filled case. 



Releasing the hard-core constraint and discussing the 
occurrence of the population imbalance in more realistic 
situations are interesting future directions. We expect 
that the basic features of the imbalanced phase and the 
transition uncovered in the present work are robust, ir- 
respective of the microscopic details. 
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